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Abstract. It is demonstrated that non-equilibrium vibrational effects are enhanced 
in molecular devices for which the effective potential for vibrations is sensitive to the 
charge state of the device. We calculate the electron tunneling current through a 
molecule accounting for the two simplest qualitative effects of the charging on the 
nuclear potential for vibrational motion: a shift (change in the equilibrium position) 
and a distortion (change in the vibrational frequency). The distortion has two 
important effects: firstly, it breaks the symmetry between the excitation spectra of 
the two charge states. This gives rise to new transport effects which map out changes 
in the current-induced non-equilibrium vibrational distribution with increasing bias 
voltage. Secondly, the distortion modifies the Franck-Condon factors for electron 
tunneling. Together with the spectral asymmetry this gives rise to pronounced nuclear 
wave function interference effects on the electron transport. For instance nuclear- 
parity forbidden transitions lead to differential conductance anti-resonances, which 
are stronger than those due to allowed transitions. For special distortion and shift 
combinations a coherent suppression of transport beyond a bias voltage threshold is 
possible. 
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1. Introduction 

The question how quantized vibrational modes affect the electron transport through a 
single molecule has recently attracted a lot of interest, both experimentally [TJ 12 E3 HI 
13 El E! and theoretically [HI EI COl EH CSl CHI CHI CHI EZl EH- For one, charging 
a molecule in a junction may induce a shuttling of the center of mass [HJ Ej through 
the interaction with the electrodes. It may also change the internal configuration of 
the molecule by changing bond lengths or angles ^Hj. In most of the cited works 
Coulomb charging effects on the molecule played an important role. Previously we 
demonstrated that this is even more so when multiple orbitals couple asymmetrically 
to a vibrational mode and start competing [TH], cf. also [2D]- Up to now theoretical 
models for transport have assumed the frequency of the oscillation to be independent of 
the charge state of the molecule. In this paper we investigate the effects of a distortion 
in addition to a shift of the nuclear potential with respect to some coordinate internal 
to the molecule. We focus on the simple case of a single orbital (i.e. two relevant charge 
states). Situations may occur where the shift of the potential is small for symmetry 
reasons and the distortion plays a prominent role. For instance for a hindered rotation, 
charging of the molecule may lower the rotational barrier without shifting the potential 
minimum. The resulting change in frequency may be large as cases are known where 
the barrier is even inverted (i.e. potential minima and maxima interchange roles |15j). 
In general, an extreme flattening of the potential surface may also lower the energy of 
a dissociative continuum of states of the nuclear motion. This will not be considered 
here, see however [HI]- We study the simplest transport model possible incorporating 
distortion effects for bound nuclear motion together with the usual shift of the nuclear 
potential minima. We emphasize that this is the most general harmonic approximation 
to the nuclear potentials of a pair of charge states involved in transport. This model 
is readily extended to more detailed potential surfaces adapted to specific situations. 
We focus mainly on low energy excitations for which anharmonic effects lead only to 
quantitative corrections. 

An important property of harmonic potentials is their spatial inversion symmetry with 
respect to the minimum. If the shift is sufficiently small relative to the distortion, the 
nuclear wave function parity is quasi-conserved. This imposes a quasi-selection rule on 
the Franck-Condon (FC) factors which determine the rates for the electron tunneling. 
This leads to effects in the intensity of the non-linear conductance resonances which 
are related to destructive interference of the nuclear wave functions. This is to be 
contrasted to the blockade [T2j HHj and NDC effects (THj discussed previously, which 
basically follow from classical features of the nuclear motion even though the discrete 
vibrational excitations are of quantum signature. Quantum effects in the intensities 
of the conductance resonances may also occur for moderate shifts of the distorted 
potentials. The FC-factors of a single low-lying vibrational excitation may become 
strongly asymmetric due to constructive and destructive interference, thereby strongly 
increasing its lifetime. In a large region of applied voltages this leads to a bias driven 
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Figure 1. Schematic effective nuciear potential (cf. Eq. J5J) for the neutral (N = 0, 
black) and charged (N = 1, blue) electronic state of the molecule. The large charging 
energy has been subtracted since near the charge degeneracy point only changes in the 
vibrational energy are of importance for the differential conductance. The potentials are 
plotted as function of the normalized coordinate Q (see text) and are relatively shifted 
by \/2X > and distorted by a factor a 2 — ujq/uji > 1. The vertical arrows indicate 
the change in the classical elastic energy when charging (blue) and de-charging (black) 
the molecule, respectively, while maintaining the nuclei at the equilibrium position of 
the initial charge state. 



population inversion of the vibrational states which suppresses the electronic transport. 
This is only possible when the zero-point motions (ZPM) of the vibration in different 
charge states are sufficiently different due to a distortion. 

The paper is organized as follows. In Sect. El we formulate the transport model and 
discuss how the distortion affects the non-equilibrium transport through the quasi- 
classical and quantum features of the FC-factors (Sect. 12. Relaxation of the 
vibrational excitations is incorporated (Sect. I2.2|) in order to identify non-equilibrium 
effects by comparison. In Sect. El we discuss the results and transport mechanisms for 
weak and strong distortion. We conclude with a discussion in Sect. |U 



2. Model 



We consider a molecule which is weakly tunnel-coupled to a left (L), and right (R) 
metallic electrode and capacitively coupled to a gate electrode. The transport is assumed 
to be dominated by a single orbital on the molecule and double occupancy of this orbitals 
prohibited by strong Coulomb blockade. For simplicity we disregard effects of electron 
spin degeneracy (the minor quantitative corrections are trivial to include) so we can 
consider spinless non-interacting electrons. Nuclear motion internal to the molecule is 
taken into account through a single vibrational mode which couples both linearly and 
quadratically to the charge on the molecule. As mentioned above, quadratic coupling 
may certainly compete with linear effects. The Hamiltonian H = Hu + Yl r ^ r + 
incorporates the molecule (M), reservoirs r = L,R and the tunneling between molecule 
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and electrodes (T): 

H M = \N)(N\H N , (1) 



N=0,1 
r? M 

Hn = ^ + y u n (* ± Sx/2f , (2) 

H e = / ] £fcrQfc r Qfcr; (3) 

r=R,L k 

r=i?,L fc 

where ± corresponds to |iV) denoting the neutral (A/ - = 0) and charged (N = 1) 
electronic state of the molecule respectively. We employ units h = ks = e = 1. The 
electrodes r = L,R are described by the non-interacting Hamiltonian H e . They are 
maintained at fixed temperature T and electrochemical potentials \L r — fi ± V/2, where 
V is the applied bias voltage. The average chemical potential is defined such that fi = 
corresponds to the charge degeneracy point (including the change in vibrational zero- 
point energy of the molecule). The gate voltage then effectively varies /i relative to 
the vibrational transition energies of the molecule. The Hamiltonians H^, N = 0, 1 
describe the vibrational motion in the neutral and charged electronic state, respectively, 
which gives rise to the excitation spectra. This is depicted in Fig. Here x and 
P ([ x iP] = "0 are the vibrational coordinate and momentum, respectively, and M is 
the effective mass. Typically, in each charge state the lowest order approximation to 
effective nuclear potential around its minimum is a harmonic potential. Importantly, 
both the equilibrium position and frequency depend on the charge state, ujq ^ u)\ and 
5x 0. Which effect will be more important depends on microscopic details of the 
molecule. Since the ZPMs are different it is not obvious how to define a single normalized 
coordinate. It is convenient to introduce the geometric-mean frequency Q = ^uj^uji 
as the vibrational energy scale. We relate the shift to the dimensionless parameter 
A = a/ MVl/25x and the frequency distortion to a = ^Jujq/uj\ (i.e. c<Jo,i = o^fi). Using 
the dimensionless coordinate Q = V MQx normalized to the ZPM associated with Q and 
the conjugate momentum P, [Q, P] = i we rewrite Hn in the form (see also Appendix A| ) 



TT U N 



P 2 + /„ . , ^\ 2 ' 
-x + a 
a* 



(5) 



with ± for A^ = 0,1. We label the molecular eigenstates by N m and write their wave 
functions as \N)\m)jf, where |m)jy is the vibrational eigenstate with m = 0,1,2,... 
quanta excited on the potential surface for electronic state \N),N = 0,1: Hjsr\m)isi = 
uj^im + l/2)|m)jy. Different regimes are characterized by comparing the change in the 
classical elastic energies involved in the vertical transitions to each potential (see Fig. H} 
with the vibrational frequency of the final charge state. (Equivalently one compares 
the ZPM of each potential with the relative shift of the two potentials.) The relevant 
dimensionless couplings are thus X 2 a ^ 1 and A 2 /a ^ 1. Due to the spatial inversion 
symmetry of each potential about its minimum the transport problem is invariant under 
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A — > —A and (a,fi) — > (a 1 , —[/,) (see Appendix A ) . Therefore, assuming A > 0,a > 1 



only three regimes need to be considered (i) A 2 < 1/ot, (ii) l/a < A 2 < a, and 
(iii) a < A 2 . 

Transport. The tunneling of electrons with an excess energy provided by the bias voltage 
drives the molecule out of electronic and vibrational equilibrium. In the weak tunneling 
regime, i.e. r < T, the occupation probabilities for A" electrons on the molecule 
and m vibrational quanta excited may be described by a stationary master equation. 
Neglecting for now relaxation of the vibrational states by the environment (see Sect. 12.2(1 
we have: 

P° = D = (w r P 1 — W r P°) 

± m w / j \ ' ' 0m<— lm' ± m! ' v lm 1 VOm 1 m ) ) 

rm! 

Pm 1 — — (Wi m '^o m P m ~~ W 0m ^ lm ,P m ,) , (6) 

rm 

with transition rates due to tunneling to/from electrode r = L,R 

^ltn'f-Om = r Pm' mf r{p-m' m) i 

W^ m ^ lm , = F r F m , m [l - fr(Hm'm)]- ( 7 ) 

where f r (E) = ^ E ~^ r ^ T + l) -1 . The addition energies for the transition l m i <— m are 
l^m'm = wim' — u;o m — ^(a -1 ™,' — am). (8) 

and the Frack-Condon factors 

F m > m = |i(m'|m) | 2 . (9) 

Where possible we will reserve m and m' for vibrational numbers of the charge state 
N = and N = 1, respectively. The stationary current flowing out of reservoir r = L,R 
is given by 

mm' 

The probabilities are normalized, YliNm ^m = 1; an( i the current is conserved, Il + Ir = 
0. 

Spectral features. Due to the /x and V dependence of the transition rates (jZJ) the current 
will change whenever a resonance fi r = fi m >m is crossed. For r = L,R this defines a line 
with negative /positive slope in (//, V) plane for > (which we consider from hereon) 
where a new electron tunneling process becomes possible. Here the molecule can change 
its vibrational energy by an amount /i m > m . Without distortion (a = 1) this resonance 
condition only depends on the change in vibrational number wl — m. For instance, once 
the transition Oo — > 1& is energetically allowed for some fixed integer k, transitions 
m — ► lm+fc are allowed for all m = 1,2,.... This infinity of processes becomes 
energetically allowed at a single resonance line. In combination with other allowed 
transitions a cascade of transitions gives access to arbitrarily high excited states jTH] 
and results in a divergence of the average phonon-number when we let A — >• [16J for 
fixed applied voltages. For a ^ 1 the mean of the two vibrational numbers also enters 
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into the resonance condition since /j, m i m /Q = (a + a -1 )^' — m)/2 — (a — a~ 1 )(m + m')/2 
The cascades of transitions are now switched on in a sequence of steps. 
Intensities. The sign and intensity of the change in the current at the resonances is 
determined by the rates T r for tunneling to/from electrode r = L,R and the Franck- 
Condon (FC) factors F m r m . The FC-factors take into account that the nuclear potential 
is altered when the molecule becomes charged. Their energy dependence through the 
vibrational numbers m, ml is typically dominant over that of the rates which we take to 
be constants r r = 2%\t r \ 2 p r with density of states p r in electrode r = L,R. In contrast to 
most transport models, here the FC-factors are non-symmetric F m / m ^ F mm i . In general 
this is the case when the nuclear potentials in the two charge states are not identical 
up to a shift. The sum rules ^2 m F m > m = J2 m , F m i m = 1 are due to the completeness 
of each vibrational basis set {|m)o} and {|m')i}. These guarantee that the current 
and the total occupation P N = Ylim °f eacn charge state N will saturate at large 
bias voltage to the electronic limit (i.e. without the vibration) I e = (^ r l/r r )~ 1 and 
P e 0,1 = F L ' R /(J2 r l/r r ), provided the FC-factors are bias voltage independent (cf. [§]). 



2.1. Franck- Condon factors - Classical and quantum features 

The FC-factors for any A and a can be calculated analytically by disentangling the 
unitary transformation which maps the oscillator Hi onto the oscillator Hq. The 
expressions and their derivation are deferred to Appendix A and we will focus here 



on their essential features. In Fig. 121 we plot the FC-factors in gray-scale for three 
representative cases. Their large-scale dependence on the vibrational numbers m, ml 
follows from quasi-classical arguments as is discussed in detail in |Appendix B| In the 
case of only a large shift (A 2 > 1, a — 1, Fig. the FC-factors are basically nonzero 
only in a classically allowed region bounded by the so-called Condon parabola [2*T| H§]. 
The maximal values occurs at m = X 2 ,m' = and m = 0,m' = A 2 . For A — > 
this parabola narrows down to a line, F m i m —> 5 m i m . For only a large distortion 
(a 3> 1,A = 0, Fig. Efc)), the classically allowed region has two linear boundaries 
(which coincide for a — > 1) and the maximal value occurs for m = m! = 0. When 
both a shift and a distortion occur (Fig. Efb)) the boundary curve is partially linear and 
partially parabolic. Additionally, two classically allowed regions of different intensity can 
be distinguished, corresponding to a difference in possible classical motions. The global 
maximum occurs for m = 0, m! — X 2 /a and a local maximum for m = X 2 a,m' = 0, 
corresponding to the two different classical elastic energy scales. It is important to note 
that for a > 1 one can have a shift which is large relative to the ZPM of one potential 
(A 2 a > 1) but still small relative to the ZPM of the other (A 2 /a > 1). In this case 
interference effects may occur in the FC-factors which are directly observable in the 
electronic transport (see Sect. 13. 2 j) . 

The distortion breaks the symmetry between the neutral and charged state in two 
respects: the vibrational excitation spectra of the neutral and charged state become 
asymmetric and the FC-factors become asymmetric as function of the energies E m = 
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Figure 2. Franck- Condon factors in grayscale (black = l,white=0). The quasi- 
classical boundary curves (yellow lines) separate classically allowed and forbidden 
transitions (see |Appendix B| . Sketched below each panel are the classical phase space 
(Q,P) trajectories of the vibration in charge state N — (black) and N = 1 (blue) 
with energy mu>o and m'u)\. The red numbers (2) and (4) in the grayscale plots indicate 
the number of intersections of these trajectories, (a) a = 1, A = 3.0. For strong shift, 
X 2 ^S> a the FC-factors are suppressed for small m',m and exponentially increase as 
the classically allowed region is approached, (b) a = 2.05, A = 3.0. For combined shift 
and distortion the FC-factors are asymmetric in m, ml (although they are symmetric 
in the limiting cases (a) and (c) due to special symmetries), (c) a = 2.05, A = 0. 
For strong distortion, A 2 -C 1/a the FC-factors rapidly decrease with m, ml and in 
addition oscillate due to the parity selection rule (checkerboard pattern). 



muJo,E m i = m'ui (not shown). Note that this is also the case for A = 0, a > 1, even 
though in this case F m > m = F mm i (see Fig. Etc)). Together these lead to a lack of 
inversion symmetry of the current with respect to the gate energy /i — > — \i even in the 
case r L = T R . 

The qualitative classical understanding of the FC-factors allows one to understand the 
dependence of the occupations and the current on the applied voltages in detail using 
figures like Fig. |21 This is discussed in Appendix C| For instance, for strong shifts 
(A 2 3> 1, a — 1) one can explain the blockade of the current at low voltage ^H], negative 
(NDC) instead of positive differential conductance (PDC), and even sharp current 
peaks ^Hj in terms of a feedback mechanism in the vibration-assisted transitions between 
the molecular states. The detailed variations of the FC-factors within the classically 
allowed region are not needed for this. Although the current steps at discrete energies 
are due to the quantized nuclear motion the variation of their sign and intensity follow 
from quasi-classical features of the FC-factors. However, when a distortion is present 
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interference effects in the sign and intensity of the conductance resonances may appear. 
For instance, for purely distorted potentials (A = 0) the spatial inversion symmetry of 
the vibrational wave function cannot be changed by the electron tunneling. This leads 
to a strict parity-selection rule for the FC-factors (see | Appendix A ): 

= unless m! — m = even (11) 

This is visible as a checkerboard pattern in Fig. Efc). For weakly shifted but distorted 
potentials the tunneling rates can still vary significantly when the vibrational number 
changes by only one (see Sect. 13. 2. If . This leads to even-odd effects in the intensity and 
sign of the conductance. For a intermediate shift of the potentials the decay rate of a 
single low-lying excited state can be coherently suppressed while its rate of population is 
coherently enhanced. This is possible only when the nuclear ZPM of the two potentials 
are sufficiently different. As a result the transport is suppressed in a large regime of 
applied voltages. Both effects require asymmetric excitation spectra i.e. a distortion. 
This indicates why interference effects in the FC-factors played no role in previous works. 

2.2. Relaxation 

The FC-factors strongly influence the type of non-equilibrium vibrational distribution 
which the transport current induces on the molecule ||. To be able to identify such 
effects qualitatively we compare our results with those where vibrational excitations on 
the molecule can relax by coupling to an environment of oscillators. For simplicity 
we assume that the spectral function of the environment is a constant 7 and the 
temperature is equal to that of the electrodes, T. For weak coupling to this 
environment, 7 <C min {u^} , its influence can be included through an additional term 
^2 n (WNnn-NnPn — WNn^NmPm) to the right-hand side of equation (jHJ) for without 
altering the expression for the current (jlOj) . The rates are given by 

W Nn ^ Nm = j[±b(u N (n - m))} (12) 

for n ^ m and N — 0, 1 where b(E) = [e E l T — l) -1 = — (1 +b(—E)). The relaxation can 
be either strong (7 > T) or weak (7 < T) relative to the tunneling as long as both are 
smaller than T and min {u>n}. What is of interest here is that non-equilibrium effects of 
different physical origin disappear at different characteristic strengths of the relaxation 
7 < T i.e. they have a specific sensitivity to relaxation processes. The case of strong 
relaxation will not be discussed except for the important fact that in this limit NDC 
effects vanish in any single orbital model regardless of the FC-factors 19J. This may 
be shown using an equilibrium ansatz for the vibrational distribution ^UJ ^] . Thus 
in limit of weak tunneling and weak coupling to the environment considered here NDC 
implies non- equilibrium. 

I We point out that for A — the model makes little sense physically without relaxation since for 
a =/= 1 the even and odd m states can not be mixed by electron tunneling processes and for a = 1 the 
vibrational number cannot change since F m > m = 5 m 'm- The master equation @ in these case has no 
unique stationary solution (i.e. the solution depends on initial conditions). This artifact immediately 
disappears when introducing a finite A or relaxation. 
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3. Results 

The stationary current / (Eq. p()jl ) and differential conductance dl/dV are presented 
for symmetric tunneling rates P* = T,r = L,R and temperature T/Q = 0.01. Gray- 
scale plots of dl/dV have different linear scale factors for dl/dV ^ to clarify voltage 
conditions for which NDC effects occur. Their magnitudes can be appreciated from the 
presented /(V) curves or from the text. We will consider parameter values which are 
well separated to keep the discussion simple. At this point a general conclusion can 
already be made: for all values of A for which we present results no NDC is visible if we 
set a — 1. The occurrence of NDC in all cases where a ^ 1 indicates that a distortion 
enhances non- equilibrium vibrational effects, however via several different mechanism 
which we will now analyze. 

3.1. Nearly symmetric excitation spectra: a 2 < 2 

For moderate values of a 2 < 2 only one vibrational excitation of AT = 1 lies below the 
first one for N = (cf. Fig. Q). The low energy spectra in the two charge states may 
thus be characterized as nearly identical. Interestingly, due to the slight asymmetry 
the current provides direct information on the changes in vibrational distribution which 
remain hidden without a distortion. Some effects of interference in the FC-factors may 
be identified. 

3.1.1. Small shift X 2 <C l/a,a: Broadening of the vibrational distribution For a = 1 
the dl/dV is completely featureless apart from the ground-state transition lines. A 
distortion causes the ground-state resonance line to split in many excitation lines which 
can be resolved at low temperature, as can be seen in Fig. El The current, shown 
in Fig. |U is strongly modulated on the new small energy scale u> — uj\ by which the 
resonances are separated. This modulation is expected since the tunneling rates strongly 
decrease with increasing energy. The featureless result for a = 1 is rather special since 
it is due to the exact symmetry of the excitation spectra for N = and N = 1. 
The progression of lines maps out the stepwise lengthening of cascades of transitions 
(Fig. Efa)). Once [l/(a 2 — 1)] + 1 ([x] = integer remainder of x) of such resonances 
have been traversed (more than 2 for a < y/2) , the cascade is infinitely long i.e. any 
state can be reached. This is the case for — //r > u . The broadening of the vibrational 
distribution is thus sharply controlled by the bias voltage in this region. The individual 
processes which lead to the broad non-equilibrium vibrational distribution for A <C 1, 
discussed in [Ug\, may thus be identified in the transport current if the vibrational 
excitation spectrum is charge dependent due to a distortion. From Fig. |2^c) one sees 
that for finite a the FC-factors do not collapse onto a line F m i m = b~ m i m as we let A — > 0. 
We therefore argue that the scaling of the vibrational distribution width found for a = 1 
in ^Bj may break down for a^l below some cut-off value for A (which depends on a). 
For 2u>i > hl > uji the progression seems to break down after the first large step. This 
is due to the suppression of all rates between even and odd excitations cf. Eq. (JTTJ), as 
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Figure 3. Differential conductance for a — 1.05, A = 0.01. White and black relative 
to the grey background indicate positive and negative intensity of dl/dV. The same 
holds for all subsequent grayscale plots. 
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Figure 4. I(V) normalized to the electronic limit I e for /i/f2 = 1.5 in Fig. 01 The 
current steps are spaced by the frequency difference luq — w\. 

Fig. E^b) illustrates. Except for one sharp feature and a weak NDC effect the resonance 
lines corresponding to quasi- forbidden transitions are missing. Pronounced NDC occurs 
when the region fi L > uji,u — u)\ > —fiR is reached: the change in the current is 1/3 
of its maximal value. The occupation of state li via a quasi-forbidden transition causes 
the drop in current since the probability is redistributed over the three states , lo 
and li and the latter does not contribute to the current since Fi F o- Once the 
transition 0i <— 1\ is energetically allowed all excited states are de-populated due to 
the quasi-selection rule for the FC- factors: F w F n . The current nearly reaches its 
maximal value T/2. For A 2 < 1/at the selection rule is sufficiently weakened that the 
progression already starts at /i^ = uj\ as seen in Fig. O 

The progression is a non-equilibrium effect as it consists only of transitions between 
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(a) Cascade switches (b) Ground-states 
on in steps carry large current 

/J L >2to 1 (o 1 </j l <2co 1 




Even excitations Odd excitations 

get populated stay depleted 



Figure 5. Transitions relevant to the right in Fig.|3 In all this kind below black/blue 
levels correspond to ground- and excited vibrational states of the neutral/charged 
molecule, TO ,m = 0,1,2,.../ l m >,m' = 0,1,2,.... The parity of the states is 
indicated by the line style (full=even, dashed=odd). Arrow thickness indicates the 
relative strength of the FC-factors. Parity forbidden/allowed transitions are indicated 
by a full/dashed arrow. Not all transitions are indicated, (a) For [i^ > 2ui% all 
grey transitions are energetically allowed and only the missing links (green) in the 
cascade need to be switched on, one by one starting from below, at the resonance 
lines — jitR = k(u!o — ui) for k = 1, 2, . . . , [l/(a 2 — 1)] + 1. After the last transition is 
possible the scheme repeats for the higher lying states i.e. all states are energetically 
accessible and a broad vibrational distribution is established, (b) For 2u>i > fiL > U\ 
the progression does not occur due to the parity effect. 
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Figure 6. Differential conductance for a = 1.05, A = 0.1. 
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Figure 8. I(V) normalized to the electronic limit I e for fi/Cl = 0.5 in Fig. Wc 
note that the average vibrational number shows the same marked peak as function of 
the bias voltage. 



excited states. The number of visible lines is indicative of the average vibrational number 
and is reduced when a relaxation rate 7 > is switched on (not shown). We point out 
that in Fig. 0] the energy separation in this progression of excitations could be mistaken 
for a vibrational frequency. Also one might infer erroneously a coupling A < 1 since 
many resonances with decreasing intensity can be resolved. In fact the extent of the 
progression shows an dependence on the shift opposite to that of a usual FC-progression: 
it becomes shorter with increasing A. Crucial for a correct identification are furthermore 
the featureless low bias region (\fi r \ < ^o^i) an d the NDC effects. 

3.1.2. Large shift l/a,a <C A 2 : Trapping in the vibrational ground state For strong 
shift the parity of the nuclear wave functions plays no role: in Fig. [7| resonance lines 
lii — mui and —fiR = mui for m = odd are now clearly visible. For a = 1 a strong 
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(a) Ground-state 
carries no current 



(b) Feedback 
not active yet 




High excitations 
decay fast 



Low excitations 
decay slowly 



Figure 9. Transitions relevant to the right in Fig.0 As in Fig. El the missing links 
(green and red) in the cascade are switched on. However, in this case the molecule 
becomes trapped in the ground states (not contributing to the current) for most applied 
voltages: once excited states become slightly occupied they feed back strongly into the 
ground state (black arrows) since the transition rates exponentially increase with the 
difference in vibrational energy (opposite to Fig. [SJ, see Fig. [21 Due to the gate 
asymmetry, fi > 0, the ground state can only decay via processes with much smaller 
rates. Due to the distortion of the spectrum this feedback mechanism is not yet active 
at the first two resonances (green arrows) , causing the current to "overshoot" resulting 
in the marked peak in Fig. [H] 



shift A 2 ^> 1 of the potentials is well-known to suppress the ground state transition 
line |22| IT2*| 113] and redistribute the weight into a FC-progression of conductance 
resonances. At this point it should be emphasized that the suppression is a non- 
equilibrium effect even though the vibrational distribution has its main weight in 
the ground states. This is evidenced by the disappearance of the NDC and current 
peaks [12] which occur for values of A > 4 and asymmetric gate energy > uj) and 
the enhancement of the suppression ^H] upon full vibrational equilibration. 
The main effect of the distortion comes from the asymmetric spectrum. What is 
remarkable in Fig. compared to Fig. |3j is that most of the resonances at higher bias 
with spacing on the small scale uuq — uj\ lead to NDC, whereas the excitations on the 
larger scale uj ,ui all correspond to PDC. The same cascade of transitions discussed 
in Sect. IH.l.ll now stabilizes the lowest vibrational states jTH] which contribute little to 
the current. However, due to the small mismatch of the excitation spectra the current 
can first reach a high value which is subsequently reduced to the value it has for a = 1 
when the cascades are switched on, see Fig. |H] This happens repeatedly with increasing 
bias. The effect here should thus not be characterized as current suppression. Rather, 
the current is enhanced relative to the case a = 1 where the feedback always dominates 
the current. Due to the slightly asymmetric spectra the trapping in the ground state is 
"postponed". Markedly, the NDC lines have positive/negative slope for ji ^ which is 
opposite to that of the NDC occurring for very large A and a = 1 |T3J EI] • 
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Figure 10. Differential conductance for a = 2.05, A = 0.01. 
3.2. Asymmetric spectra: a 2 > 2 

For sufficient distortion, ujq > 2u\, two or more excitations of the N = 1 charge state 
lie below the first excitation in the N = state i.e. there is a true asymmetry between 
the spectra of the two charge states. This brings about a simplification: most of the 
features discussed below are qualitatively reproduced when truncating the spectrum at 
energies above the larger vibrational frequency ujq, retaining only the states m — 0, 1 
and m! = 0, . . . , [a 2 ] + 1 for N = and N = 1 respectively. Due to the presence of 
the N — 1 low-lying states interference effects in the FC-factors also gain importance. 
For small shift the quasi-conservation of the nuclear wave function parity suppresses the 
electron tunneling between all even and odd vibrational states of iV = and N — 1. 
For larger shifts A 2 > 1 /a interference may suppress the decay rate of a single low-lying 
excitation and simultaneously enhance the rate at which this state is populated. This 
concerted effect of constructive and destructive interference is due to the opposite parity 
of the N = vibrational ground- and excited state. 

3.2.1. Small shift X 2 < I /a: Parity effect The dl/dV in Fig. QUI has two distinctive 
features. Below the threshold voltage where 0i can not yet be reached (horizontal zig- 
zag line), all resonance lines (with negative slope) correspond to NDC with an intensity 
which alternates. Above the threshold resonances lines (with positive slope) appear 
with alternating sign of the differential conductance. Both effects derive from the parity 
quasi-selection rule incorporated in the FC-factors and are systematic: as one increases 
a 2 by one, the additional resonances appearing below and above the threshold voltage 
follow the above pattern. We explain this parity effect using Fig. [T2l 
NDC intensity alternation. All the low voltage resonances correspond to NDC since 
the low- lying excitations of N — 1 become successively occupied (Fig. fT2T a)) without 
contributing significantly to the current (the FC-factors F m i rapidly decrease with ml). 
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Figure 11. I(V) normalized to the electronic limit I e for fi/Cl — 0.1 in Fig. The 
first three current plateaus values are describe by Eq. 1(13)1. In particular, for a 2 ^> 1 
the initial current drop by a factor 1/3 signals the redistribution of the probability 
from two to three states. 



(a) NDC intensity 
alternation 



(b) Conductance sign alternation 
PDC NDC 






Odd excitations 
contribute less 



£?_ 

Odd excitations Even excitations 
feed back into even feed back into odd 



Figure 12. Transitions relevant to the parity effects in Fig. ^] The closely spaced 
excitations belong to charged state N = 1, the levels with large spacing (only two 
shown) belong to N — 0. (a) Successive occupation of the non-contributing states l m ' 
leads to NDC lines at low bias. The variation of the NDC intensity reflects wether the 
transition is parity forbidden (dashed line) or allowed (full line), (b) The successive 
opening of escape processes (green and red arrows) from the states not contributing 
to the current gives PDC for even and NDC for odd excitations. 



Since the rates at which each excitation is populated and depleted are the same (although 
small), the ground state occupation is reduced oc 1/V resulting in a lower current at 
higher bias. However, the FC-factors F m 'o also oscillate: transitions from the ground 
state to odd-m' states l m i are strongly suppressed relative to those with even m', 
F{2n+i)o F(2 n )Q for n — 1, 2, — Remarkably, the quasi-forbidden transitions appear 
as anti-resonances in the differential conductance (instead of missing resonances) which 
modulate the current stronger than allowed transitions. This effect is related to the 
non-equilibrium conditions (see below). This may be explicitly demonstrated from the 
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Figure 13. Differential conductance for a — 2.05, A = 0.01 and finite relaxation 
rate 7/r = 0.1. At this intermediate value of the relaxation only the odd NDC lines 
are suppressed, cf. Fig. The spacing between the remaining even resonance lines 
with negative slope may be mistaken for the vibrational frequency of the excitations i.e. 
they mimic an effective frequency doubling. However, the resonance lines with positive 
slope are inconsistent with such an interpretation and hint at the missing resonances 
due to the parity effect. 

expression for current plateau k — 1, 2, . . . where k oc V 

k 

h = [(k + i)/r R + i/r L ] - 1 ^ F m , . (13) 

m'=0 

This result applies at low bias where states m = and states m' < k are occupied and 
state 0i is not yet reachable (Fig. H2fa)). For T r = T the current after an initial big step 
decreases oc 1/V, since the number of occupied excited levels which do not contribute 
to the current grows oc V. Interestingly, for asymmetric coupling to the electrodes 
r L <C T R the depletion of the ground states is postponed: the current initially increases 
in steps with V to reach a maximum around V oc uj\{T R /T L ) and then decays oc 1/V. 
For negative bias the current shows no such maximum. 

NDC / PDC alternation. From the low bias resonances one can directly find the FC- 
factors if the tunneling rates T L,R are known. However, above the threshold voltage this 
is not the case anymore: here multiple states from both charge sectors contribute in a 
more complicated way (Fig. IT^T b)). Now escape from the non-contributing states l m > 
becomes possible. Remarkably, this suppresses the current further at lines with positive 
slope terminating at the NDC resonances due to quasi- forbidden transitions. This is 
due to a feedback mechanism in the vibration assisted transitions which effectively traps 
the system in the odd-parity states as explained in Fig. El This is somewhat similar to 
the mechanism in the opposite case of strong shifts ( Sections 13.1.21 and I3.2.3J) but relies 
critically on the modulation of the rates due to the quasi-conservation of parity. 
The allowed and forbidden excitation lines have different sensitivity to relaxation 
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processes and disappear in two stages. When increasing the relaxation rate 7, in a first 
stage the strong NDC effects due to quasi-forbidden transitions become comparable in 
intensity with the NDC due to allowed ones and subsequently disappear as shown in 
Fig. [T31 Thus similar to optical spectroscopy the parity selection rule now leads to 
missing resonances in the spectrum. The alternation of the sign of the resonances above 
the threshold voltage also disappears since it is caused by asymmetries in the smallest 
rates. Only in a second stage the remaining NDC lines due to the mf =even states 
disappear. 

3.2.2. Intermediate shift 1/a < A 2 < a: Coherent suppression Interestingly for any 
a > \/2 a drastic suppression of the current occurs near special values of A, a prominent 
example being 



As seen in Fig. H31 at large voltage bias the current is completely suppressed beyond 
the ground-state transition line /jlr = 0. In this region the low-lying excited state I2 
becomes completely occupied i.e. we have a bias driven population inversion between 
vibrational states (cf. JH])- The reason for this is twofold and is explained in Fig. 115( b): 
due to destructive interference the rate of decay to the ground state 0o is suppressed 
(FC-factor F20 = for A = A(2)(a;)); simultaneously constructive interference maximizes 
the tunneling rate into this state from the excited state 0i. This concerted effect is 
due to the opposite parity of the ground- and excited state for N = 0. As soon as 
the excited state 0i can be reached via some tunneling processes the excited state 
1 2 becomes fully occupied and suppresses the transport. This happens at the bias 
voltage threshold forming the horizontal zig-zag line which we encountered above. The 
transport is recovered only when direct escape (li <— O2) from the coherently blocked 
state becomes energetically allowed i.e. —fiR > uj — 2uji (strong white line with positive 
above the suppressed region in Fig. [TJJ. The effect is coherent in the sense that both 
destructive and constructive interference of the nuclear wave function are responsible 
for the suppression of electron transport through the molecule. 

In a similar way, the FC-factor F m > of a higher excited state l m /,m' = 3,4, •• • may 
vanish for some value of A. (For mf — 1 this happens only for the trivial value A = 0.) 
If in addition this state lies below the first excitation for N — 1, i.e. a 2 > mf, this leads 
to a region of suppressed current similar in shape to that in Fig. El but more narrow 
(e.g. for mf — 3 the width of the region is halved). The lines in the (a, A 2 )-plane where 
both these conditions for the coherent suppression are satisfied are plotted in Fig.lTolfor 
mf = 2 — 9. The curves are all of the form Eq. (fHj) but with a prefactor which differs 
from 1/2. The region where this interference effect occurs is centered around the regime 
1/a < A 2 < a and a > ^/2, where it is possible to have a shift which is larger than the 
ZPM of the flattest potential but still smaller than that of the steepest potential. With 
increasing a the values of A where the coherent suppression occurs start to abound and 
even proliferate into the regime where the shift becomes strong, A 2 > a, see Fig. [H3 




(14) 
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Figure 14. Differential conductance for a = 2.05, A = 0.7. 




Figure 15. (a) Transitions relevant to the coherent suppression of the current on the 
right in Fig.^J The decay rate of excited state ±2 to the ground state 0o is suppressed 
due to destructive interference in the nuclear wave function overlap. State I2 can thus 
not be reached directly, but due to the asymmetric spectra it can reached in three 
tunneling process via the excitations above it and the excited state 0i (for V above 
the zig-zag threshold in Fig. I14f> . The rate for the last process, 0i — > I2 is maximized 
due to constructive interference. The coherently enhanced ratio of rates for going in 
and out of I2 suppresses the current, (b) Nuclear wavefunction overlap for A = A( 2 ): 
the maximum and node of the states Oo and lo, respectively, align with the node of 
excitation ±2, thereby creating a low-energy blocking state. 

With increasing a the number of such values rapidly increases roughly oc a 4 : since state 
l m > has m! nodes there are [m'/2] zeros of F m /Q as a function of A for a given a and only 
the states m! — 2, • • • , [a 2 ] have energy < uj$. Finally we note that excited states with 
m' > [a 2 ] are not expected to cause a coherent suppression effect since they always have 
two (groups of) states with opposite parity to decay to. It is highly unlikely that the 
decay rates to both types of states can be suppressed simultaneously by a special choice 
of the shift A. 

Naturally the coherent effect is more sensitive to parameter values than the quasi- 
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Figure 16. Conditions for coherent suppression in the (a, A 2 ) plane. Regimes of weak 
(A 2 < 1/a), intermediate (1/a < A 2 < a), and strong shift (a < A 2 ) are separated 
by black boundary lines. The blue curves labeled by integers m' indicate the values 
of A for which the FC-factor F m /Q vanishes for fixed a. The curves are plotted only 
for those values of a for which state m! lies below the first excitation of N = charge 
state and where the coherent suppression occurs. 

classical trapping effect (cf. Sect. 13 . 1 7^1 and ET.2.rSj) . This sensitivity has an interesting 
side to it. (The effect of voltage dependence of A on the trapping effect was considered 
in jSj.) Introducing only a weak dependence of the parameters, for instance A, on the bias 
voltage, the current exhibits a pronounced dip down to zero when tuning the parameters 
a and A with V through the condition for coherent suppression (fTlj) . I n view of the many 
situations where this effect can occur (Fig. ITH|) this is an interesting novel possibility to 
be explored in single molecule devices. 

Similar to the parity effect, introducing relaxation affects the transport in two stages. 
First the suppressed groundstate transition line is restored and the forbidden transition 
line disappears as depicted in Fig. H3 An NDC effect related to the state with the 
suppressed FC-factor is still present. In a second stage this effect also vanishes. Excited 
states m! < [a 2 ] have effective decay rate ~ (wf — 1)7 and are therefore more sensitive 
to relaxation. In summary: of all the excitations l m > the state m! = 2 gives rise to 
the strongest coherent suppression effect in the largest voltage region and is the least 
sensitive to relaxation. 

3.2.3. Large shift A 2 > a: Trapping in the vibrational ground state The strong 
asymmetry between the excitation spectra due to the distortion results in the 
pronounced asymmetric conductance plot in Fig. As one approaches the strong shift 
regime the NDC at resonances hl = m'uj\,m' = 1, . . . , [a 2 ] discussed below (Fig. fTUj) 
turns into PDC (cf. the first excitation in Fig. H3J) and becomes suppressed in intensity 
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Figure 17. Differential conductance for a = 2.05, A = 0.7 and relaxation rate 
7/T = 0.1. Note that the excitation line of state with suppressed FC-factor is now 
missing and that the NDC at the zig-zag is weakened and moved to the straight line 
/jl = 3wi. 
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Figure 18. Differential conductance for a — 2.05, A = 3.0. Note that the first 
excitation now leads to PDC, cf.Fig.0 

at low bias for A 2 > a as seen in Fig. EH This is also described by Eq. (|13|) (which 
holds for any A) since the FC-factors increase with m',m for A 2 > a, cf. Fig. E[b). 
For negative gate energy n the current slowly increases, but for positive /i the current 
surprisingly shows a sharp increase followed by NDC, despite the strong shift. Fig. I27H 
explains how the postponement of the classical trapping effect (similar to that discussed 
in Sec. I3.1.2J) allows for large current steps despite the strongly shifted potentials. This 
results in strong PDC excitations spaced by the larger frequency ujq and in between 
several pronounced NDC excitations with smaller spacing uj\ (in contrast to Fig. [7| 
where the NDC spacing is ojq — ui). 
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Figure 19. I(V) normalized to the electronic limit I e for fi/fl = 1.25 in Fig. 1181 
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Figure 20. Transitions relevant to the NDC effects on the right in Fig. (a) 
For most bias voltages transitions are switched on (red arrows) which eventually feed 
back (black arrows) into the vibrational ground state lo • Escape from the highest 
states [a 2 ], [a 2 ] — 1, ... with largest contributions to current becomes possible first, 
leading to the strongest NDC. (b) At the ground state transition line all the excitations 
l m i,m' — 1, . . . , [a 2 ] initially become occupied (gray arrows) and contribute to the a 
relatively large current without the feedback occuring. 

4. Discussion 

We have found that non-equilibrium vibrational effects are enhanced in molecular 
devices for which the effective potential for vibrations is sensitive to the charge state 
of the device. We modeled this by a change in the vibrational frequency in addition 
to a shift of the potential minima. In particular, for weak distortion of the potential 
the current was shown to map out sharp changes in the vibrational distributions with 
bias voltage. For sufficiently strong distortion of the potential interference effects of 
the nuclear wave functions were show to strongly influence the electron transport. The 
coherent effects are also expected to occur for more detailed models of the nuclear 
potentials since the requirements on the low-energy vibrational excitations are rather 
basic. The parity effect requires potentials which are distorted and only slightly 
shifted. The coherent suppression due to the vanishing overlap between an vibrational 
excited- and ground-state corresponding to a different charge on the molecule, requires a 
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moderate shift and distortion. The precise conditions for the suppression will be different 
but Fig. El provides the qualitative picture. These mechanism, together with possible 
weak dependence on current or voltages offer interesting possibilities for controlling 
electron transport in single molecule devices. 



Appendix A. Franck-Condon factors for shifted and distorted potentials 



In this Appendix we give the expressions for the FC-factors for potentials exhibiting both 
a relative shift and distortion. The derivation can be done by straightforward algebra 
without recourse to special functions (cf. [2H1 12H ) • We first note that the sign of the 
shift is irrelevant since the each of the states in the overlap integral has a definite parity 
with respect to spatial inversion relative to the minimum of its potential. Using this 
we find that interchanging u <-» u)\ (i.e. a = ^JZj^Juj[ — ► ft -1 is equivalent to charge 
conjugation N ^ 1 — N, N = 0, 1 or // — > —fi. We can thus restrict ourselves to a > 1 
as long as we discuss both polarities of \x. For the calculation of F m i m it is convenient 
to normalize the coordinate to ZPM of the potential in question, Qn = x\JMuj^ and 
its conjugate 



H 



N 



UN 

2 



P N + (Qjv ± Va~±/2\ 



(A.1) 



where N = 0, 1 corresponds to ± and the shift parameter is A = a/ 'MVl/25x and the 
mean frequency Q = y/cJoUh- We obtain oscillator 1 from oscillator by first applying 
a shift and then a distortion: 
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We write the corresponding transformation of the lowering operators, 
(Q N + iP N ) /V2, 
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unitary transformation b\ = e s boe s , where 
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The exponential can be disentangled by the methods described in 
e L6t 2 +zfet e (inC)(b t 6+2)+c e _Rfe 2 +rfe w h ere parameters in the exponential factors are 
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For a = 1 one obtains the well-known expression F^r^ = e x2 ™^\ 2 \ m m '\ 
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where L is the associated Laguerre-polynomial and m < m! and F n 
A = the selection rule Eq. (fTT]) is easily verified. The nonzero matrix elements, for 
instance, for the special case m = 0,m' = 2k are F 2 ko = yl" 



e { ^e k where 

(2ife)!! = {2k){2k - 2) • • - 2 , (2fc - 1)!! = (2k - l)(2Jfe - 3) • • • 1 and ^ = ttFt- The 
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expression in terms of the Hermite polynomials H n reduces to the known result for 
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Appendix B. Classical features of the FC-factors 

The large scale variations of the FC-factor F m ' m in the (m, m') plane and their effects on 
the transport have a simple classical interpretation. It is important to discuss these if 
one wants to identify quantum effects of the nuclear motion. The central point is that the 
FC-factor becomes exponentially suppressed unless the nuclear motions in the effective 
potentials of the two charge states are compatible i.e. the phase-space trajectories of the 
two motions intersect. The boundary between classically forbidden and allowed regions 
in the (m, m')-plane is found by requiring that the simultaneous equations (cf. Eq. (JSJ)) 

aP 2 + -(q + ^) =2m' (B.l) 



a V y/2, 

have at least one real valued solution for (P, Q) (the vibrational energies are Eq = lo^m 
and E\ = oj\m'). Within the classically allowed region there may be regions of different 
overall intensity related to the appearance of additional solutions. The intersections of 
the elliptic orbits determined by Eq. (|B.1|) are illustrated in Fig.0 In the case of shifted 
potentials, A > 0, a = 1, Q always has a single real solution. Two real solutions for P 
exist if 



m + m > - 
~ 2 



\2 , 1 / /\2 

A + - m-m 
A z 



(B.2) 



This parabola in the (m, m!) plane, tilted by an angle 7r/4 relative to m! axis, is the 
so-called Condon-parabola [21J depicted in Fig. Ufa) . The condition (|B.2|) is equivalent 
to demanding that the classical turning points of the two motions are interspersed. In 
the case of distorted potentials, A = 0, a > 1, the requirement is 

1 in' n t v 

-a < — < c? B.3 
cr m 
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and real solutions are always four in number. The left inequality ensures that Q is 
real. The corresponding lower boundary line in Fig. WL C ) is equivalent to requiring the 
potential energies of the motions are equal at the maximal coordinate (P = 0). The 
right inequality ensures a real solution for P and corresponds to equal kinetic energy at 
Q = (i.e. not at the classical turning point). For the general case A > 0,a > 1 two 
real solution exist for P when m, m! lie inside a Condon-parabola 



which is tilted by the angle tp relative to m axis, where tany? = a~ 2 = uji/ujq i-e. the 
parabola tilts towards the axis corresponding to the highest frequency. It touches the 
axes at m = X 2 a and m' = X 2 /a, respectively, corresponding to the elastic energies. The 
solutions for Q are always real in this region. However, an additional two real solutions 
for Q occur between the parabolic boundary and below the line 



beyond its tangent point with the parabola, which is located at m — (A 2 a)/(a 4 — 
l) 2 ,m' = (A 2 /a)/(l — a~ 4 ) 2 . In Fig. Efb) one sees that in this region the FC-factors are 
clearly enhanced. The above classical expressions thus give an simple guide to the large 
scale structure of the complicated exact expression (jA.2|) . For A — > the parabola (jB.4|) 
becomes very narrow and reduces to a line through the origin line m'/m = 1 /a 2 , i.e. 
we recover Eq. (|B.3|) . For a — > 1 the region with 4 solutions moves to infinity with the 
tangent point of the parabola and we retain (jB.2|) . 

Appendix C. Vibrational distribution and numerical convergence 

Much of the behavior of the occupations P^ and the current as function of the applied 
voltages may be understood from a simple scheme which is readily extended to situations 
with multiple competing orbitals [E2]. First, to determine which direct transitions 
between states relevant one draws the region llr/VL < m'/a — ma < lil/Q ("bias 
window") into the grayscale plot of the FC-factors. For m',m within this region 
transitions l m > <-> m are both allowed whereas above/below this region only l m > — > m 
I l m i <— m is allowed by electrons leaving / entering the molecule via both tunnel 
junctions. Importantly, m', m outside the classically allowed region can be disregarded. 
In a second step we have to determine which states will become occupied significantly 
via cascades of tunnel processes (unless the relaxation is extremely fast, 7 3> T, in which 
case we can solve Eq. (JBJ) using the vibrational equilibrium ansatz J01E])- These may in 
principle allow arbitrarily high states to be reached. The FC-factors which satisfy sum 
rules (cf. EI) will prevent the average vibrational numbers ^m m -^ > m fr° m increasing 
indefinitely with V. The expressions for the boundary curves (jB.4|) and (jB.5|) now 
become helpful in estimating the number of vibrational states required to solve the 
transport problem with good accuracy. The points of intersection (m'.,m r ) of the edges 
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of the bias window with the boundary curves can be explicitly found (e.g. for a — 1 
these points have a quadratic dependence on both V and We may truncate the 
infinite set of master equations beyond the cut-offs for on m and m! estimated as 

m c = max r {m r (/i, V)} , m' c = max r {m' r (/i, V)} . (C.l) 

Beyond these points the asymmetry between the FC-factors in the gain and loss terms 
in the stationary master equation increases with m and m! . Therefore the occupations 
will start to strongly decrease. The convergence of the current requires more states 
to be taken into account for strong shifts A 2 a ± ^> 1 even when the distribution is 
already converged. In this case the exponential increase of the FC-factors and the 
strong decrease of the occupations with m tend to cancel out. For small shift the cut- 
offs m c , m' c overestimate the number of required states. The distortion generally widens 
the classically allowed region i.e. transitions with larger change m! — m become more 

probable which improves the convergence. 
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